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Abstract. We study a large class of models with an arbitrary (finite) number 
of degrees of freedom, described by Hamiltonians which are polynomial in bosonic 
creation and annihilation operators, and including as particular cases n-th harmonic 
generation and photon cascades. For each model, we construct a complete set of 
commuting integrals of motion of the Hamiltonian, fully characterize the common 
eigenspaces of the integrals of motion, and show that the action of the Hamiltonian 
in these common eigenspaces can be represented by a quasi-exactly solvable reduced 
Hamiltonian, whose expression in terms of the usual generators of 5(2 is computed 
explicitly. 
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1. Introduction 

Dating back at least to 1990, when Zaslavskii |1J pointed out the relation between some 
quasi-exactly solvable (QES) systems 0, ||, 0, EL [7| and the Dicke and Heisenberg 
models, there has been an ongoing interest in the correspondence between the exactly 
calculable part of the spectrum of QES systems and some finite-dimensional systems 
like the spin systems discussed by Zaslavskii or the effective Hamiltonians frequently 
used in nonlinear optics. 

In essence, Zaslavskii's method consisted in making a suitable ansatz for the 
coefficients of the eigenvectors of each "multiplet" of the spin system, coefficients which 
could be grouped into a generating function that turned out to be related by elementary 
transformations to the wavefunction of the QES system. 

In 1995, Alvarez and Alvarez-Estrada || studied the usual second harmonic 
generation effective Hamiltonian by showing its equivalence to an infinite family of one- 
dimensional QES systems. In their method, the transition from the "discrete" photon 
system to the "continuous" QES systems is achieved by first transforming the effective 
Hamiltonian from the second quantization to the Bargmann representation. The key 
idea to separate variables and identify the one-dimensional QES systems — the second 
step of their procedure — is to use as a variable an appropriate quotient of powers of 
the Bargmann variables that describe each oscillator. (Later, the same authors || used 
this method to study the problem of third harmonic generation and pointed out that it 
works in general for the n-th harmonic generation.) 

Two recent papers by Dolya and Zaslavskii [|TI], [DJ proceed in a related but different 
direction: they study models with a single degree of freedom whose Hamiltonians are 
even polynomials in the creation and annihilation operators, and show that under certain 
conditions these Hamiltonians also lead to a QES system. 

Finally, an important idea for our work can be traced back to the early papers 
on QES models by Ushveridze and Zaslavskii, and appears with more or less emphasis 
in [|], H ||, [□], |[: the role played by integrals of motion to reduce completely integrable 



Hamiltonians with more than one degree of freedom to a family of QES systems. 

In this paper we study a large class of models with an arbitrary (finite) number of 
degrees of freedom described by Hamiltonians which are polynomial in the creation 
and annihilation operators. These models include as particular cases the effective 



Hamiltonians of n-th harmonic generation [12, 13], as well as photon cascades of equal 



or different frequencies ||14|| , which, to the best of the authors' knowledge, have not 
been analyzed before as QES systems. For each model, we first construct a complete 
set of commuting integrals of motion of the Hamiltonian. We use the Bargmann 
representation to completely characterize the common eigenspaces of the integrals of 
motion, which are of course invariant under the Hamiltonian. We show that the action 
of the Hamiltonian in these common eigenspaces can be represented by a QES reduced 
Hamiltonian, whose expression in terms of the usual generators of is computed 
explicitly. We emphasize that the derivation of this explicit expression is not just an 
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exercise of academic interest but the starting point for the application of asymptotic 
methods to calculate the corresponding eigenvalues in the limit of a large number of 
photons. By way of example, we give explicitly the expressions for n-th harmonic 
generation and iV-photon cascades, which in the former case can be compared with the 
particular instances n = 2, 3 studied in || [|. The paper ends with a brief summary 



2. The models 



In this paper we shall study the following general Hamiltonian 

H = H + gH 1 (1) 



where 



N M 



H = J2 u i a l a i + J2^ b l b k V,Hk>0 N,MeN (2) 
i=i k=i 

M N N M 

Hi = n( & D mfc n a i + u^ ni ■ n ^ ^ ^ e n ( 3 ) 

k=l 1=1 1=1 k=l 

and the frequencies v\ and are subject to the energy-conservation constraint 

N M 

^nivi = y^m fc /i fc . (4) 
1=1 k=\ 

In the preceding expressions, the operators ai,b k represent boson annihilation 

(creation) operators with frequencies vi,fj, k , respectively. The Hamiltonian ([]]) is an 
effective model describing the conversion of a number of photons of frequencies u±, . . . , 
into photons of frequencies fit, ... , \im- It includes as particular cases the usual processes 
of n-th harmonic generation (if N = M = mi = 1, m — n) and multiple photon cascades 
(if M = mi = ni = ■ ■ ■ = njv = 1), which have recently received considerable attention 



in the literature [T2|, [13], 14" | 



A key property common to all these models is that they are completely integrable in 
the sense that there exists a set of iV + M (the number of degrees of freedom) pairwise 
commuting, functionally independent operators including the Hamiltonian. These sets 
are clearly not unique. We have found it convenient to take as our set 

{H ,H u A h B k | / = 1,...,JV-1; k = 1, . . . , M - 1} (5) 

which includes the unperturbed Hamiltonian H and the perturbation Hi separately, 
and where 

Ai = m + ia\ai - nia] +1 ai+i (6) 

B k = m k+ ib[b k - m k b\ +1 b k+ i. (7) 

A straightforward computation using the standard bosonic commutation relations shows 
that these operators commute pairwise, and we shall see below that they are in 
fact functionally independent. As a consequence, there exists a basis of common 
eigenfunctions of all the operators in the set. We will see in the next section that 
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it is possible to give an explicit description of the simultaneous eigenspaces of H Q , A t 
and Bk (the "unperturbed eigenspaces"). By studying the action of Hi within these 
subspaces, in the next section we shall establish the exact solvability of the model (|l])-(|4]) 
and derive some general properties of its spectrum. 



3. Invariant subspaces and matrix representation 

Following H, we shall study the spectrum of the Hamiltonian (|l|) using the Bargmann 
representation [15]. In this representation, the Hilbert space for a single boson is the 
space of entire functions of the form 

oo 

n vnl 

n=0 

where the complex numbers c n are such that 

oo 

^|c„| 2 <oo (9) 



n=0 



and the scalar product is defined by 

1 



7T 



(g, f) = - d(Rez)d(lmz)g(z)f(z)e-^ . (10) 



The orthonormal harmonic oscillator eigenstates \n) are given in this representation by 

z n 

|n)--= n = 0,l,... (11) 
Vnl 



and the annihilation and creation operators are represented by derivation with respect 
to z and multiplication by z, respectively 

d + , s 

a -> — a) -)• z. 12 
dz 

The Hilbert space TC of the model ([]]) is the tensor product of iV + M single boson 

Hilbert spaces. Denoting by xi (I = 1, . . . , N) and y& (k = 1, . . . , M) the complex 
variables associated to the "a" and "6" degrees of freedom, respectively, we have the 
following assignments: 

ai -> d Xl a] ->• xi b k -> d Vk b\ -> y fc . (13) 

A basis of the Hilbert space Ji of our models is the set of (in general unnormalized) 
monomials 

xY i e Zf, j G (14) 

where Z + denotes the set of nonnegative integers, and we have used the multiindex 
notation 

N M 

x^Hxf y j = ]Jyi k - (15) 

1=1 k=l 
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The corresponding expressions for the operators in the set ([5]) are 

N M 

Ho = 1 ^2 v ^ x ^i + VkWedyk ( 16 ) 

1=1 k=l 

H l = y m d n x +x n d™ (17) 
A t = n l+1 x t d Xl - mxi +1 d Xl+1 (18) 
B k = m k+1 y k d yk - m k y k+1 d yk+1 (19) 
where x n and y m are monomials as in equation ( jl5|) and 

N M 

% = ^ = U d 2- ( 2 °) 

Z=l fc=l 

(In the rest of the paper we shall make frequent use of the multiindex notation without 
further notice whenever there are no ambiguities in the interpretation of the formulae.) 

We first note that each monomial ([14]) is a common eigenfunction of the operators 
H , Ai, B k , I = 1, . . . , N — 1, k = 1, . . . , M — 1, with eigenvalues respectively given by 

N M 

E o = ^2 Uiii + ^ k ( 21 ) 

1=1 k=l 

ai = n l+1 ii - n t i l+1 l = l,...,N-l (22) 
(3 k = m k+1 j k - m k j k+1 k = I, ... ,M — I. (23) 

The spectrum of the Hamiltonian (|IJ) can thus be computed by diagonalizing the 
perturbation Hamiltonian H 1 in the common eigenspaces SE ,ai,...,a N . 1 ,/3 1 ,...,/3 M . 1 , spanned 
by the monomials (H) whose exponents satisfy equations (f2~TD-([23|) for fixed values of 
Eq, oti, . . . , Pm-v 

We shall now describe the common eigenspaces 5^ iai) ... lQ!JV ._ 1)( g 1) ... )( g M _ 1 = S. Since 
the unperturbed energy Eq is a finite sum of nonnegative terms, it follows that S is 
finite-dimensional. Thus, the model ([]]) is exactly solvable, since its whole spectrum can 
be computed algebraically. 

Let x p y q be a given (fixed) monomial in S. We shall first prove that S is spanned 
by monomials of the form 

fs = xPy% s (24) 

where 

and s ranges over a finite interval so < s < s\ of Z. 

Indeed, let x p y q be any other monomial in S. The quotient Q = 

X P -Pyl ~q i s 

then 

a generalized monomial (i.e. a monomial with possibly negative integer exponents) and 
thus may not belong to TC. However, the action of the operators (|Ttj|)-(|T9"|) on such 
generalized monomials is well-defined, and a simple calculation shows that Q satisfies 
the system of first-order partial differential equations 

HqQ = AxQ = ■■■ = A N -iQ = BxQ = - - • = B M -iQ = . (26) 
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In other words, Q must be a joint invariant of the M + N — 1 commuting vector fields 
H , Ai, B k . A second straightforward computation shows that the function ( defined 
in (|25|) is a joint invariant of this set of vector fields. Since the number of independent 
variables is N + M, the general solution of the system ( p6|) is an arbitrary (smooth) 
function of (, and since Q is a generalized monomial, it must be a power of £, which 
proves the first part of our claim. 

Since S is finite-dimensional, the values of s for which f s lies in S must be bounded 
above and below. If 

So = min{s G Z : f s G S} S\ = max{s G Z : f s G S} (27) 

we shall next prove that 

S = (f s | s < s < si) (28) 
where ( ) denotes the linear span. Indeed, since clearly 

H f s = E f s A l f a = a l f a B k f s = (3 k f s (29) 

for all integer values of s, it suffices to verify that for all s G [so, si] fl Z the exponents 
of f s are all nonnegative integers. Since s G [s , s x ] and both / so and f 81 belong to S, 
this follows from the inequalities 

Pi + ms >pi + n t s > q k — m k s > q k - m k si > . (30) 

This completes the proof of our claim. 

Note that (|2~7D and the inequalities ( |3~0D imply that the bounds s and Si are given 

by 

s = m 



(- 


El 




. n i. 



{ 


' Qk ' 




_m k _ 



mm 



: 1 < k < M } (34) 



: 1 < / < X} .si ^minj-^- : 1 < k < Mj (31) 

where [■] denotes the integer part. Calling 

M = Pi + s ni > M k = q k - s m k > (32) 
the eigenspace S can be alternatively written as 

5e5J = I V(1,(,...,0 (33) 

where 

Mk 
m k _ 

and the nonnegative integers Mi and M. k are subject to the single restriction 

Mi < ni for al least one I G {1, . . . , iV} (35) 

which is an immediate consequence of the definition of sq. The set of N + M nonnegative 
integers {Mi, ■ ■ ■ ,Mn,MIi, ■ ■ ■ ,-Mm} subject to the condition ( |3~5"|) define uniquely the 
eigenspace S. Indeed, Mi is the minimum power of x\ and M. k is the maximum power 
of y k of the monomials in S. We have thus proved the main result of this section: 

Theorem 1 The common eigenspaces of the operators H , Ai, B k , I = 1, . . . , iV — 1, 



k = 1, . . . , M — 1, are the spaces given in equations (j33|)-(|3 
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Since these eigenspaces Sj^ are invariant under H 1 (because this operator commutes with 
H and Ai, B k for all I, k), we devote the next paragraphs to study the corresponding 
action. 

Consider the basis B of Sj^ spanned by the normalized vectors 

1 x x +snyM - sm s = ,l,...,r (36) 



\[C L 

where (see equation ([□])) 

N M 

C s = + n lS )\ ■ H(M k - m k s)\ . (37) 

i=i k=i 

A straightforward computation using equations (|i"7|), ([3~6"1) and (|37| ) yields 

fr — N 



s i=i 



M 



+ J ^ - m ^ s ) • • • (A^fc - m fc (s + 1) + 1) e s+1 . (38) 

fc=i 

Thus, the matrix Hx representing i^x in the basis B is tridiagonal, Hermitian, and has 
zero diagonal entries. The only nonzero entries of Hi are given by 

i 

(N ni—1 M m k -l \ 2 

J] + ^8 + ^ + 1) fj (M k - m k s - i k )\ =(H 1 ) a>8+1 (39) 

J=l J;=0 fc=l i fe =0 / 

with s = 0, 1, . . . , r — 1. 

A number of immediate conclusions can be drawn from the structure of the matrix 
Hi. We note in the first place that Hi acts irreducibly in <Sj^. This implies that Hi 
is functionally independent from H , A\ and B k for all k, I. Since x n y m is annihilated 
by all the A\ and B k but not by Ho, the latter operator is functionally independent 
from the former ones. Since all the Ai and B k are clearly functionally independent, the 
preceding arguments show the functional independence of all the operators in the set fl5|) 
and finish the proof of the complete integrability claimed at the end of section 2. 

Secondly, the eigenvalues of Hi are real, simple, and symmetrically distributed 
around zero. Indeed, let S S (E), 1 < s < r + 1, be the s-th principal minor of the matrix 
E~K — Hi. A straightforward computation shows that 5 S (E) satisfies the three-term 
recursion relation 

5 s+ i (E) =E5 S (E) - h 2 s 5 s ^(E) s > 1 (40) 

where h s = (Hi) SjS _i and 5i(E) = E, 8 (E) = 1. Since h 2 s > for 1 < s < r, lemma 1 
(section 1.8) of [16] implies that the roots of the characteristic polynomial 5 r+ i(E) are 
real and simple. They are also symmetrically distributed around zero, for 5 r+ i(E) con- 
tains only either even or odd powers of E on account of the form of the recursion 
relation (M). 
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Remark. The model ([I]) admits certain solvable (and in some cases integrable) 
generalizations. In the first place, we can replace gHi by any polynomial in Hi with 
real coefficients. Clearly, the resulting Hamiltonian is completely integrable, since it 
commutes with the operators (|5|). Moreover, this generalized Hamiltonian preserves the 
finite-dimensional subspaces Sj^ and is therefore exactly solvable. 

More generally, one can consider a perturbation Hamiltonian H[ consisting of a 
finite sum of terms of the form (§), provided that the exponents in each term satisfy 
condition (gj). In other words, 

(M N N M \ 

n^-n^+n^-nc (41) 
k=i 1=1 1=1 k=i j 

where S is any finite subset of the set of all pairs (n, m) e Z% x z + satisfying the 
condition (|J). Note that Hq still commutes with H[, since it commutes separately 
with each term in the sum fl4l[) . Thus the spectrum of H can be computed by 
diagonalizing the new perturbation Hamiltonian H[ in the eigenspaces of Hq, which 
are finite-dimensional by equation (^l|). Hence this more general model is also exactly 
solvable. Note, however, that now to each term in the sum ([|l|) there correspond two 
families of operators A\ = Ai(n) and Bk = Bk(m) which do not necessarily commute 
with the remaining terms in H[. Hence, this model need not be completely integrable 
for general values of the coupling constants g nm . 



4. Quasi-exact solvability of the reduced Hamiltonians 

In the previous section we have shown that the computation of the spectrum of the 
Hamiltonian ([I]) is equivalent to the diagonalization of a family of tridiagonal Hermitian 
matrices with zero diagonal elements. For small values of r + 1 = diniiS^, the 
eigenvalues of the matrix Hi representing Hi can be easily computed in closed form. In 
many situations of physical interest, however, the number of photons of each frequency 



involved can actually be very large. This implies (see equation (p4|) ) that the parameter 



r can also take large values. Although the spectrum of Hi can be determined numerically 
for fixed values of the parameters Mu M-k labelling the space Sj^, one is often interested 
in the behavior of its eigenvalues as a function of A/], M-k (and thus r). In the case of 
second-harmonic generation, i.e. N = M = mi = l,ni = 2, Alvarez and Alvarez- 
Estrada p| derived asymptotic formulae for the eigenvalues of Hi by performing a 
semiclassical analysis of the second-order ordinary differential equation (ODE) obtained 
by restricting the perturbation Hamiltonian Hi to SjJl = S*, where e = 0, 1. It was 
also shown in that paper that the restriction of Hi to S 1 ^ is in fact equivalent to a QES 
operator in a single variable, in agreement with previous results due to Zaslavskii |T|. 

In this section we show that the restriction of our general perturbation Hi to an 
eigenspace Sj^ is also equivalent to a QES differential operator in a single variable. 
Besides the purely mathematical interest of this result, the knowledge of explicit 
expressions for the ordinary differential operator of the restriction Hi\ s j^ is precisely 
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the starting point to apply asymptotic techniques, which are currently being developed 
for third- and higher-order ODE's. 

Consider the action of H\ on an element x Jsf y M P(( > ) G Sj^ 

M N N M 

k=l 1=1 1=1 k=l 

where P is a polynomial of degree at most r. Since 

n i — i 

W*MO) = ^ _n ' II CM - ii + C»c)^(C) (43) 

mfe-l 

^ fc (rf^(0) = 2/f fe - m& II - - m * ^C)^(C) (44) 

*k=0 

it follows that 

(xV^(O) = x^y M H 1>Ted P(C) (45) 
where the reduced Hamiltonian is given by 

V n;— 1 M m k -l 

H ltied = - ii n cm - ji + ^ c«9c) + c n n (-^^ - ^ - mfc w (46) 

^ i=l j i= fc=l i k =0 

Thus the action of H\\ s u is equivalent to that of Hx jZed on the space V r of polynomials 
in £ of degree at most r. Note in particular that the invariance of <Sj^ under implies 
that V r is invariant under Hi ied (this fact can be also verified directly from equation ( PH) 
using ( |34]) and (^); see the discussion below). 

Let d = max(^ 1 n h Y^k=i m k) be the order of the differential operator H l red . If 
d < r, a well-known theorem due to Turbiner [^7| (see [IB for a simplified proof) states 
that the invariance of the polynomial space V r under H l ied implies that this operator 
can be written as a polynomial of degree d in the generators of the realization of 
spanned by 

J + = ( 2 d ( -r( J^ = Cd ( - r - J~ = d c . (47) 

If d > r, Turbiner's theorem only guarantees that the r-th order part of Hi iied is 
a polynomial of degree d in the operators fl4"T|). In fact, even if d > r the reduced 
Hamiltonian Hi )md can be written as a polynomial of degree d in the generators (|47D . 
Indeed, choose k' G {1,...,M} and I 1 G {1,...,N} such that (see equations ([34]) 
and (H)) 

my 

Since all the factors in the products of the expression ( f46"D commute, it follows that 

— n f*^*i)- n n>+^+0 

^=0,5,,^,, V 7 i=l,l^i'i,=0 V 7 



and A/]' < n v . (4£ 
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my — l 

+ (-m) m J + - n ( J 

i k i=0,i k i^A4 k i mod my 
M m k -l 



h' -My r 

171k' 2 



■ n n^+^+i)- w 

k=l,k^k' i k =0 V 

This general expression of the differential operator Hi^ ed as a polynomial in the 
generators of the 5^ algebra is of course not unique, and can be written in different forms 
using the commutation relations of the generators (f4T|). We conclude this section by 
showing explicitly the physically most important particular cases, in which the general 
expression fl49|) simplifies considerably. 

Example 1 Consider the problem of n-th harmonic generation, in which N = M = 
mi — 1, ni = n (and therefore the multiindices reduce to ordinary indices). The 
eigenspaces of H are S^. = x e y r V((), where r e Z+, e = 0, . . . ,n — 1, and ( = x n /y. 
The corresponding expression of Hi Tcd reads 

H hved = n n J~ H (j° + e —l + r -j-J+. (50) 

This expression can be easily shown to be equivalent to the expressions obtained in || 
for second harmonic generation (n = 2) and in for third-harmonic generation (n = 3) 
and e = (incidentally, in the cases e = 1 and 2 were not considered). 
Example 2 In the case of a multiple photon cascade, M = mi — n\ — ■ ■ ■ — njsr — 1, 
the eigenspaces of H and A;, I = 1, . . . , JV — 1, are = x-^y r V((), where My = for 
some V G {1, . . . , iV}, r G Z + , and C = (%i ■ 4 -x^/y. The expression of Hi red in terms 
of the generators (0) is 

N 

Hi, cA = J~ J] (j°+A^ + 0-J+. (51) 
i=i,i^i' 



5. Summary 

In this paper we have established a common framework to deal with a large class of 
processes in nonlinear optics including the problems of n-th harmonic generation and 
multiple photon cascades. The distinguishing feature of these models from the physical 
point of view is the condition of energy conservation, which is the key mathematical 
condition to prove their complete integrability and the exact solvability. By using the 
Bargmann representation, we have been able to provide an explicit description of the 
common unperturbed eigenspaces, as well as equally explicit expressions of the action of 
the perturbation Hamiltonian in these unperturbed eigenspaces. We have also derived 
some general properties of the spectrum of the restriction of the perturbation Hi to each 
unperturbed eigenspace, such as the nondegeneracy and symmetric distribution of the 
perturbation energies around zero. The link between these finite matrix representations 
and the corresponding continuous QES system is most conveniently established by 



QES models in nonlinear optics 



11 



the introduction of a new "projective" coordinate, which is a quotient of powers of 
the Bargmann variables that describe the different oscillators (and physically carries 
information on the phase difference among the oscillators). Furthermore, we have been 
able to give explicit expressions of the reduced QES Hamiltonians as polynomials in the 
generators of the standard QES realization of sfe by first-order differential operators. 
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